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Summary 

A  contact  fatigue  life  analysis  was  performed  for 
a  constant  speed  ratio,  Multiroller  Planetary 
Traction  Drive.  The  methodology  was  based  on  the 
Lundberg-Palmgren  theory  for  rolling-clement 
bearing  life.  The  effect  of  stress,  stressed  volume, 
and  depth  to  the  critical  stress  were  taken  into 
account.  A  design  traction  coefficient  of  0.05  was 
maintained.  Life  adjustment  factors  due  to  advances 
in  rolling-cletnetit  bearing  materials,  lubrication  and 
design  were  considered  as  well  as  the  potentially 
adverse  effects  of  traction.  Elastohydrodytiamic  film 
thickness  calculations  were  made  for  representative 
operating  conditions.  The  drive  consisted  of  a  single- 
stage  planetary  configuration  with  tw'o  rows  of  five 
stepped  planet  rollers  each  between  a  concentric  ring 
roller  and  suti  roller  element.  These  rollers  formed  a 
drive  cluster  approximately  0.21  meter  in  dlimeter, 
0.06  ntetci  in  width,  weighing  9  kilograms.  Ihc  drive 
speed  ratio  was  14.7  to  1. 

The  drive  system  10-percent  life  based  oti  tn.r 
Lundberg-Palmgren  theory  modified  by  rolling- 
element  contact  fatigue  life  adjustment  factor^  was 
calculated  for  a  spectrum  of  operating  conditions 
The  system  life  ranged  from  18  800  hours  at  16.b 
kilowatts  (22.2  hp)  and  25  (XX)  rptn  sun  roller  .speed 
to  305  hours  at  maximum  load  and  speed  of  149 
kilowatts  (200  hp)  and  75  000  rpm.  At  a  nominal 
continuous  power  and  speed  rating  of  74.6  kilov.aits 
(KX)  hp)  and  75  (KK)  rpm,  the  drive  system  life  was 
2440  huuis.  The  shortest  lived  comact  was  the  first 
planet  roller-second  planet  roller  contact  since  it  had 
the  highe.st  maximum  Hertz  stress. 

A  study  of  the  effect  of  first  and  second  row  roller 
diameters  and  center  locations  was  performed  for  a 
given  de.sign  traction  coefficient  with  a  fixed  number 
of  planet  rollers  and  for  contacts  with  constant 
transverse  geometries.  For  a  given  ratio  an  optimum 
life  geometry  existed. 

From  this  analysis,  torque  capacity  was  found  to 
be  proportional  to  the  2.8  power  of  overall  drive  size, 
for  a  given  drive  configuration  at  a  constant  traction 
coefficient  and  fatigue  life. 


Introduction 

Traction  drives  have  been  of  interest  to  inventors, 
engineers,  and  designers  for  many  years  (ref.  I).  The 


essential  feature  of  traction  drives  is  the  transfer  of 
mechanical  power  using  steel  elements  in  rollitig 
contact  wiih  the  traction  forces  being  transmitted 
through  a  thin  film  of  lubricant.  The  simplest  form 
of  It  action  drive,  although  unlubricaied,  is  the 
locomotive  wheel  upon  a  steel  rail.  There  are  dozens 
of  traction  drive  designs,  distinguished  by  the 
geometric  shape  and  arrangement  of  .ne  traction 
eiemenis  (refs.  I  to  4).  Recent  adva'-.^es  in  steel  and 
lubricant  technology  have  greatly  enhanced  the 
traction  drive  concept.  Cleaner  steels  (vacuum 
induction  melt-vacuum  arc  remeli)  have  improved 
fatigue  life  almost  an  order  of  magnitude  (ref.  5). 
Lubricants  with  a  50-percent  increase  in  available 
traction  coefficient  have  been  formulated  (ref.  4). 
Roller-io-roller  loading  can  be  decreased,  resulting  in 
longer  life  for  traction  drives. 

Traction  drives  have  many  advantages.  They  are 
quiet  and  smooth  in  operation.  In  one  case  it  was 
shown  tlta  a  traction  drive  was  24  decibels  quieter 
.hail  .r  comparable  gear  drive  (ref.  6).  High  reduction 
ratio  drives  with  high  speed  inputs  are  possible.  It  is 
also  possible  to  have  traction  drives  with 
continuously  variable  ratio  adjustment.  Traction 
drives  can  be  made  competitive  with  gear  drives  in 
efficiency,  weight,  and  size. 

Commercially  available  traction  drives  are 
manufactured  in  the  United  States,  Europe,  and 
Japan.  The  drives  arc  used  mainly  for  industrial 
machinery  applications,  but  raiely  exceed  11 
kilowatts  (15  hp)  rating.  Traction  drives  have  been 
experimentally  investigated  for  automobile 
transmissions,  submersible  ocean  going  vessel 
propulsion  drives,  and  aircraft  auxiliary  power  unit 
drives.  Workers  in  the  field  are  considering  future 
applications  to  land,  marine,  and  air  transportation 
(ref.  3). 

While  there  are  detinite  advantages  with  traction 
drives,  there  are  also  areas  of  concern  to  designeis. 
Roller  stability  has  been  studied  (ref.  7). 
Maintenance  of  adequate  lubricant  film  thickness  is 
desirable  to  prevent  rapid  wear  of  the  traction 
surfaces.  Sufficient  loading  on  the  traction  elements 
is  necessary  to  pi  event  destructive  gross  slipping. 
High  contact  stress  due  to  load  on  the  rollers  limits 
the  surface  fatigue  life  of  the  drive.  Each  of  these 
factors  should  be  careful'y  addressed  when 
evaluating  a  proposed  design  for  a  traction  drive. 

Recently  a  multiroller  planetary  traction  drive  was 
tested  at  the  Lewis  Research  Center.  It  is  shown  in 
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figure  1.  The  drive  performed  smoothly  with  a  peak 
loading  of  127  kilowatts  (170  hp),  73  cioo  rpm  input 
speed  and  nominal  14.7  to  1  speed  ratio.  The  peak 
efficiency  was  95  percent  (ref.  8).  The  multiroller 
design  concept  was  presented  in  reference  9.  This 
type  of  drive  has  nigh  power  capability,  inherently 
balanced  internal  forces,  making  bearing 
requirements  less  severe,  and  offers  very  high 
reduction  ratios. 

XhP  niirnn^/*  /'tf  thic  rpnort  is  aivp  an  anal^rcJc  f/%r 

the  life  of  multiroller  traction  drives. The  analysis  is 
based  on  Lundberg-Palm.gren  theory  which  is  used  to 
predict  the  fatigue  lives  of  ball  and  roller  bearings 
(ref.  10).  Coy,  et  al  conducted  an  analysis  similar  to 
the  present  investigation  for  the  toroidal-type 
traction  drive  (ref.  11).  The  method  was  successfully 
applied  to  spur  and  helical  gears  (ref.  12)  as  well  as 
forming  the  basis  of  fatigue  life  estimations  for 
remanufaciured  aircraft  bearings  (ref.  13). 


Symbols 

a  semimajor  axis,  m  (in.) 

B  second  row  planet  roller  position 

t>  seiniminor  axis,  m  (in.) 

C  first  row  planet  roller  position 


c  orthogonal  shear  stress  exponent 

E  modulus  of  elasticity.  N/m^  (psi) 

e  Weibull  exponent 

F  curvature  difference 

/  roller  width,  m  (in.) 

h  depth-to-critical  stress  exponent,  m  (in.) 

K  constant  of  proportionality 

L  life,  hr 

/  length  of  rolling  track,  m  (in.) 

Q  normal  roller  load,  N  (lb) 

/  radius,  m  (in.) 

5  probability  of  survival 

u  stress  cycles  per  revolution 

y  stressed  volume,  m^  (in. 3) 

w  semiwidth  of  rolling  track,  m  (in.) 

z  depth  10  maximum  orthogonal  reversing 

shear  stress,  m  (in.) 

p  curvature  sum,  m  '  (in.->) 

n  life  (millions  of  stress  cycles) 

p  available  traction  coefficient  of  lubricant 

contact 

P*  design  traction  coefficient 

/'  Poisson’s  ratio 

0  maximum  contact  pressure,  GPa  (psi) 

7  maximum  orthogonal  reversing  shear  stress, 

GPa  (psi) 

'.*>  angular  velocity,  rpm 

Subscripts; 


A,B  elastic  bodies 
p  planet 

r  ring 

S’  sun 

x.j/,-:  reference  planes 

1  sun 

2  first  row  planet  (driven  track  in  speed  reducer 

mode) 

3  first  row  planet  (driving  track  in  speed  reducer 

mode) 

4  second  row  planet  (driven  track  in  speed 

reducer  mode) 

5  second  row  planet  (driving  track  in  speed 

reducer  mode) 

6  ring 
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Analysis 

Fatigue  Life  Mode! 

In  reference  10,  Lundberg  and  Palmgren  presented 
their  hypothesis  for  the  statistical  relation  of  rolling- 
element  bearing  life  to  load  and  bearing  size.The 
mode  of  failure  was  assumed  to  be  subsurface 
originated  fatigue  spalling.  The  theory  has  been 
vvide'y  used  for  bearing  life  calculations  and  also 
provides  a  basis  for  gear  life  calculation  which  has 
agreed  very  well  with  the  data  from  gear  life 
endurance  tests  (ref.  14).  The  important  parameters 
are  number  of  stress  cycles  ?),  magnitude  of  critical 
stress  7,  amount  of  stressed  volume  K  and  depth 
below  the  surface  at  which  the  critical  stress  occurs  z. 
The  stressed  volume  is  taken  as 


y  oc  wz/ 


(I) 


where  /  is  the  length  of  the  rolling  track  which  is 
traversed  during  one  revolution  of  the  rolling  body 
and  the  semiwidth  of  the  rolling  track  is  designated 
w. 

From  reference  10,  the  probability  of  survival  S, 
for  a  bearing  contact  is  given  by  the  following 
expression: 


log^oc 


rCrjfP 


(2) 


This  relation  is  consistent  with  experimental 
observations  in  the  case  of  fatigue.  The  formula 
reflects  the  known  fact  that  the  more  localized  the 
stress  is  in  the  material  (less  stressed  volume),  the 
greater  is  the  endurance.  This  is  because  on  a 
statistical  basis,  there  is  less  likelihood  of  a  fatigue 
nuclcation  site  being  coincident  with  a  condition  of 
high  stress.  Conversely,  there  is  a  greater  probability 
of  a  crack  forming  in  the  zone  of  maximum  critical 
stress,  because  the  material  is  more  rapidly  cycled 
toward  failure  in  that  reg'on.  Hence  the  depth  to  the 
critical  stress,  as  well  as  the  magnitude  of  the  stress 
important,  and  with  each  stress  cycle  the  probability 
of  failure  increases. 

From  reference  11,  the  number  of  stress  cycles 
endured  with  90  percem  reiiabiliiy  is  given  by  the 
following  equation: 


»)  = 


(3) 


Based  on  life  testing  of  air-melted  sttel  lolling- 
clement  bearings,  the  following  values  aie  valid  for 
equation  (3):  A:=1.43x10’5  (S.I.  units),  3.58x  1056 


(english  units);  =7/3;  c  =  31/3;  and  e-  10/9  (point 
contact),  3/2  (line  contact). 


Analysis  of  Single  Contact 

Figure  2  shows  two  elastic  bodies  in  roiling 
contact.  This  figure  is  a  general  representation  of  any 
traction  contact.  The  roller  crown  profiles  lie  in  the 
transverse  plane,  which  contains  the  axes  of  rotation. 
The  rolling  radii  determine  the  nominal  velocity  ratio 
of  the  two  bodies.  The  plane  of  rolling  is 
perpendicular  to  the  transverse  plane.  The  centers  of 
the  principal  radii  of  curvature  are  contained  in  the 
intersection  of  the  transverse  plane  and  the  rolling 
plane.  After  the  bodies  are  pressed  together,  the 
contact  point  is  assumed  to  flatten  into  a  small  area 
of  contact  which  is  bounded  by  an  ellipse  with  major 
axis  2a  and  minor  axis  lb.  The  pressure  distribution 
is  semiellipsoidal  as  shown  in  figure  3.  To  complete 
the  life  calculation  of  equation  (3),  the  contact 
stresses  must  be  obtained.  The  following  procedure 
may  be  used. 

A  widely  used  method  of  describing  the  geometry 
of  two  ellipsoidal  solids,  A  and  B,  in  contact  is  to 
express  it  in  terms  of  the  effective  curvatures,  that  is, 


1  -  JL  i. 

R~  Ry 

where 

JL- + 

R.\  ^Ax  '"Bx 
11^1 


(5) 

(6) 


Figure?.  -  Geometry  for  bodies  in  rolling  contact.  Conve.x  radii  are 
positive;  concave  radii  are  negative. 


3 


7 


figure  3.  -  Contact  pressure  distrihution.  Two 
general  curved  surfaces  in  point  contact. 


The  variables  /Jj^and  represent  the  effective  ladius 
of  curvature  in  the  principal  x  and  y  planes, 
respectively.  It  is  assumed  that  conve.x  surfaces 
exhibit  positive  curvatures  and  concave  surfaces, 
negative  curvatures.  Planes  x  and  y  arc  the  respective 
planes  of  maximum  and  minimum  relative  curvature 
for  the  bodies.  These  planes,  called  the  principal 
planes,  are  mutually  perpendicular.  Planes  x  and  y 
must  be  chosen  so  that  the  relative  curvature  in  plane 
X  IS  greater  than  in  plane  y,  that  is. 


’’By 


The  radii  of  curvature  may  be  positive  or  negative 
depending  on  whether  the  surfaces  are  convex  or 
concave,  respectively. 

After  the  bodies  are  pressed  together,  the  contact 
point  is  assumed  to  flatten  into  a  small  area  of 
contact  which  is  bounded  by  an  ellipse  wiih  major 
a,xi.s  2a  and  minor  axis  2b  as  shown  in  figure  3.  Plane 
”  coniains  the  major  axis  of  the  contact  ellipse  and 
plane  x  contains  the  minor  axis.  The  ratio  k  =  a/b  is 
called  the  eilipticity  ratio  of  the  contact.  The  values 
of  k  range  from  1  to  a>  for  various  curvature 
combinations  of  contracting  surfaces.  For  cylinders 
in  contact,  the  eilipticity  ratio  is  co,  and  the  flattened 
area  of  contact  is  a  rectangular  strip.  For  spheres  in 
contact,  the  eilipticity  ratio  is  I.  The  first  type  is 
called  line  contact  and  all  other  types  are  called  point 
contact. 

The  curvature  difference  is  defined  as 

Brewe  and  Hamrock  (ref.  15)  provide  simple 
formulas  for  determining  the  elliptic  integral  of  the 
first  T  and  second  E-  kind  as  well  as  the  eilipticity 
parameter  k 


5  =  1.528  +  0.602  In-/ 


f;=  1.000  + 


(Ry/R,) 


A:  ^1.0.34 


These  simple  equations  eliminate  the  need  for 
numerical  evaluation  on  a  digital  computer.  These 
simple  equations  are  then  directly  used  to  evaluate 
tne  seminiajor  a  and  semiminor  b  axis  of  the  contact 
ellipse  and  the  deformation  at  the  center  of  the 
contact  b.  Thus, 

‘'=LV2€^/V^/  J 

where 

_ 2 _ 

The  maximum  reversing  orthogonal  shear  stress 
occurs  at  a  depth  Zp  underneath  the  surface  of 
contact.  Expressions  for  these  parameters  can  be 
obtained  from  Lundberg  and  Palmgren  (ref.  10)  as 


V2f-1 

^o- 


2t(IA  1) 


(r+-l)v'27^ 


In  these  equations  t  is  an  auxiliary  parameter  and  is 
related  to  k  as  follows  (given  by  private 
communication  with  Brewe  and  Hamrock): 

/=  1  +  0.304  A-i-8?6  (17) 

In  general,  the  width  of  the  rolling  track  is  identical 
to  the  width  of  the  contact  ellipse.  However,  when 
the  width  of  the  contact  pressure  ellipse  is  exactly 
equal  to  the  width  of  the  roller,  point  contact  is,  in  a 
practical  sense,  transformed  into  line  contact.  With 
that  loading  where  point  contact  changes  into  line 
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coiuaui,  ilic  piesbuic  area  is  still  elliptical  in  shape 
and  the  press"rc  distiihtition  is  still  parabolic.  But 
with  increased  loading,  a  leveling  of  pressure  occurs. 
If  the  rollers  in  contact  are  of  unequal  width,  then 
with  increased  load  the  highest  contact  pressures  exist 
at  the  roller  ends.  This  is  the  so-called  edge  loaded 
condition.  There  is  no  known  exact  meihod  of 
calculating  the  tiue  stress  at  this  condition.  As  a 
practical  matter,  edge  loaded  conditions  should  be 
avoided. 

If  the  traction  contact  consists  of  two  rollers 
without  crown  radii,  and  having  equal  widths,  pure 
line  contact  is  said  to  exist.  Pure  line  contacts  are 
generally  avoided  in  practice  since  high  stresses  can 
occur  at  the  edges  of  the  contact  with  even  a  modest 
amount  of  misalinement.  The  line  contact  pressure 
distribution  is  shown  in  figure  4.  Tor  this  special 
case,  the  following  equations  arc  used; 


7  =  0.;s  (20) 

Z-=0.5  b  (21) 

To  account  for  the  difference  in  the  pressure 

distribution  between  point  and  line  contact,  the 
stressed  volume  should  be  taken  as  follows  (ref.  10): 

y=  \flz  (22) 
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Figure  4.  -  Contact  pressure  distribution.  Iwo  txxlles  in  line 
contact. 


Application  to  Multiroller  System 


After  the  stress  conditions  of  each  contact  location 
in  a  traction  drive  arc  determined,  the  life  in  millions 
of  stress  cycles  for  each  roller  clement  is  computed 
from  equation  (.1).  The  life  in  hours  is  given  by  the 
following  equation,  where  u  is  the  number  of 
identical  stress  cycles  per  revolution  and  w  is  the  rpm 
of  the  roller: 


Wu?  60 


(23) 


Since  the  planet  rows  do  not  orbit,  the  rotational 
speeds  of  the  rollers  arc  inversely  proportional  lo  the 
roller  diameters.  From  probability  theory,  the  life  ot 
the  system,  Z.,  consisting  of  N  roller  elements  is 
related  to  the  life  of  each  element  by  the  following 
equation: 


. (i)”- 


1  (24) 


This  form  of  the  equation  is  uselul  if  the  values  of  the 
Weibull  parameter  are  not  equal  for  each  traction 
contact,  as  would  be  the  case  for  mixtures  of  line  and 
point  contact.  If  all  exponents  are  equal,  the 
equation  is  simpler  to  .solve. 

L  =  {Lf  +  L-{+  .  .  .  +L-j0'''^  (25) 


Results  and  Discussion 

Application  to  Multiroller  1  raction  Drive 

Figure  5  is  a  schematic  of  a  two  roller -row-type 
drive  as  described  in  reference  9.  In  the  case  of  a 
speed  reducer,  the  input  torque  acts  on  the  sun  roller 
shaft  and  power  is  transmitted  through  the  first  and 
second  row  rollers  to  the  ring  roller  winch  is  the 
output  element.  Other  variations  (kinematic 
inversions)  of  the  drive  are  possible.  In  this 
publication  the  drive  was  analyzed  as  a  speed 
reducer. 

In  order  to  transmit  the  traction  forces,  the  rollers 
must  be  pressed  iogethe:  under  load.  This  is  done 
with  a  roller  ramp-type  loading  mechanism.  When 
the  drive  is  placed  under  torque  loading,  the 
mechanism  automatically  adjusts  the  roller  contact 
normal  loads  in  proportion  to  the  transmitted 
torque.  The  ratio  of  traction  force  to  normal  load,  or 
the  design  traction  coefficient  g',  is  held  constant  by 
the  loading  mechanism  for  all  torque  loads  above  a 
preset  minimum.  However,  due  to  the  static  force 
balance  that  exists  across  the  roller  cluster  as 
generated  by  the  automatic  roller  loading 
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Figure  b.  -  Cross  section  of  Nasvytis  muitiroHer  traction  drive. 


niechanisin,  the  ratio  of  traction  force  to  applied 
nonnal  load  or  design  traction  coefficient  /<*  may 
differ  on  each  roller  contact .  In  order  to  prevent  slip, 
the  design  Liaciion  coefficient  /<*  at  each  contact  is 
selected  to  be  less  than  th.e  available  traction 
coefficient  of  the  lubricated  contact  ft.  This  is 
accomplished  by  choosing  the  proper  cone  angle  in 
the  roller  loading  mechanism.  The  available  lubricant 
traction  coefficient  depends  on  many  factors,  such  as 
speed,  load,  temperature  a^  well  as  contact  geometry 
and  surface  finish,  l  or  the  purpose  of  this  analysis  it 
was  assumed  that  the  lubricant’s  available  traction 
coefficient  significantly  exceeded  the  design  traction 
coefficient  of  0.05  under  all  operating  conditions. 
Parametric  tests  (ref.  8)  performed  on  a  Nasvytis 

iilLilLiiuliti  uiiVc  OVCi  iht  OpCrSUng  COnuuU''iiS 

consideied  here  with  a  traction  fluid  established  that 
this  is  a  correct  assumpiion. 

The  Nasvytis  multirollei  drive  which  was  analyzed 
for  life  and  contact  stress  had  five  rollers  in  each  of 
two  rows.  The  test  drive  roller  cluster  was 
approximately  0.21  meter  in  overall  diameter  and 
0.06  meter  in  width,  weighing  9  kilograms.  It  had  a 
geometrical  speed  ratio  of  14.7  to  1.  Calculations 
were  made  to  a  maximum  power  level  of  149 
kilowatts  (200  hp)  at  a  maximum  sun  roller  speed  of 
75  000  rpm.  Rolling  diameter  contact  ratios  of  the 
drive  tested  in  reference  8  for  ihc  sun  roller -first 
row  roller  (driven  track)  contact,  first  row  roller 
(driving  tract) -second  row  roller  (driven  track) 
contact,  and  second  tow  roller  (driving  track) -ring 
roller  contact  were  1.28,  3.87,  and  2.97,  respectively. 

Stress  analysis  showed  that  the  maximum  contact 
pressure  or  Hertz  stress  for  the  drive  at  any  torque 


load  occurred  at  the  contact  between  the  first  ,and 
second  row  rollers.  At  a  sun  roller  torque  of  19.0 
N-m  (168  in.-lbO  this  ma.xiniuiii  Hertz  stress  was  ’.70 
GPa  (246  000  psi).  The  maximum  Hertz  stress  for 
each  contact  over  the  range  of  operating  sun  roller 
torques  is  shown  in  figure  6. 

Life  Analysis  Without  Life  Adjustment  Factors 

The  Luiidberg-Palmgren  life  analysis  presented 
previously  was  used  to  determine  the  10-perceni  (90 
percent  survival)  fatigue  life  of  each  traction  contact 
of  the  drive  tested  in  reference  8  as  a  function  of  sun 
roller  torque.  Equation  (16)  was  used  to  statistically 
sum  up  the  individual  contact  lives,  in  order  that  an 
overall  drive  system  life  could  be  computed.  Values 
of  C  appropriate  for  line  and  point  contact  were  used 
as  required.  The  results  of  this  analysis  are  shown  iti 
figure  7.  It  is  instructive  to  note  that  the  life  of  the 
total  drive  system  can  never  be  greater  than  th.e  life  of 
the  shortest  lived  contact.  It  is  also  clear  from  the 
slope  of  the  straight  line  relation  shown  on  the  log- 
log  coordinates  of  figure  7  that  fatigue  life  is 
inversely  proportional  to  the  cube  of  torque.  That  is 


The  life  given  in  figure  7  is  based  on  the  original 
Lundbeig-Palmgien  analysis,  developed  about  30 
years  ago.  The  constants  and  exponents  of  variables 
used  in  their  analysi'  were  empirically  determined 
from  life  tests  of  bearings  of  that  time  period.  Since 
that  period  there  has  been  significant  progress  in 
bearing  material  development  and  processing  as  well 
as  a  better  understanding  of  the  operating  variables 
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ihat  govern  rolling-element  fatigue  life.  For  these 
reasons  fatigue  lives  predicted  using  classical 
l.undberg-Palingreii  theory  tend  to  be  overly 
conservative  and  need  to  be  adjusted  for 
advancements  made  in  rclling-elenient  bearing 
technology.  Appropriate  life  adjustment  factors  are 
discussed  in  a  later  section  of  this  report. 

To  determine  the  effect  of  size  (roller  diameter  and 
contact  width)  on  relative  fatigue  life  at  constant 
torque,  excluding  life  adjustment  factors  and 
changes  in  traction  coefficient,  all  roller  dimensions 
were  scaled  by  a  uniform  amount  and  the 
caiculaiions  repeated.  The  results  of  this  exercise  are 
shown  in  figure  8.  It  is  clear  that  relative  life  of  the 
drive  system  is  quite  sensitive  to  size,  in  fact. 


Sue  scale 
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proportional  to  the  8.4  power  of  size.  That  is,  for  all 
elliptical  contact  shapes, 

Loc(size)**  (27) 

for  coiiSian;  torque  and  traction  coeffici-eni.  The 
value  of  the  size  exponent  will  change  if  theie  is  a 
mixture  of  line  and  elliptical  contacts  in  the  drive. 
Calculations  also  showed,  that  for  a  given  drive 
configuration,  iraciion  coefficient  and  fixed  service 
life,  torque  load  can  vary  as  a  function  of  overall 
size.  For  all  elliptical  contacts  from  equations  (26) 
and  (27),  tins  relation  is 

Toe  (size)’  ** 

as  shown  in  figure  8. 

I  rom  equation  (27),  a  14  percent  increase  in  drive 
size  triples  the  calculated  drive  life  at  the  same  torque 
level.  These  calculations  were  made  assuming  no 
change  in  applied  iraeiion  coefficient.  Flowever,  a  14 
percent  increase  in  size  at  constani  torque  and 
traction  coefficient  will  cause  a  14  percent  increase  in 
contact  velocity  and  a  33  percent  reduction  in  contact 
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picssmi;  Km  ati  cllipiically  sluipcd  coiuaci.  This  will 
uiuioiil'icdly  lead  lo  a  icdiiciioii  in  the  available 
lubiicani  tiaeiiiin  eoel  1  ieicni .  In  luin,  a 
ennipensaiinti  ineieave  in  the  applied  iiaeiion 
ei'eltieieiu  will  be  letp.iiied  Menee,  an  ineieawin  the 
ei'tttaet’s  noiinal  load  will  iwetn .  T  he  tiet  lesiilt  will 
in  tael  leduee  tlie  si/e  exponetit  in  t!ie 
aloieiiieiilioned  evination  by  bome  amount  depemlini^ 
on  the  settsiti\ity  ol  the  hibiieaiuK  iiaetinn 
eoeitieieiu  to  the  eltaiiites  iti  operating  eoiidiiions 
eited.  1  o;  example,  it  a  14  peieent  ineieabc  in  xi/c  loi 
a  given  lubricant  at  a  paiiietikn  opeiattng  speed  and 
toiqne  eaiises  a  Id  peieent  loss  in  the  available 
tiaetiiMi  eoettieieni,  drive  life  will  be  2.3  times  loiigei 
instead  ol  three  times  longet  vvheti  ebanges  in 
ttaetiiMt  eoettieieni  are  igiuncd. 

L  ite  Adjustnieni  I'aetois 

Adv  aneenient  s  in  nilling-elemcnt  bearing 
teehnology  have  prodrteed  bearing  designs  whieh,  I'oi 
a  given  applivation,  have  ineieased  fatigue  lives.  A 
better  undei standing  of  the  variables  that  affeei 
fatigue  life  and  the  use  of  itnproved  materials  and 
niatiufactui ing  iccliniqucs  bavc  icsulted  in  rolling- 
eleinetu  beaiitig  lives  iltat  ate  often  many  times  those 
whieh  Luttdbcrg-Paltngreti  theory  ptediets.  In 
rceognition  of  these  impi oveinciits  life  adjustment 
faetois  have  been  developed  in  referenee  5  for 
pioviditig  a  tnotc  aeeuratc  prediction  of  rolling- 
element  bearing  life.  In  view  of  the  similarity  in 
contact  geotnetry,  tnaterials,  lubricating  factors, 
operating  conditions  and  failure  tnodes,  several  of 
tliesc  factors  developed  feir  rolling-element  bearings 
ate  eoiisidered  to  be  equally  applicable  to  ihc 
multiroller  it  action  diive.  The  life  adjustment  factors 
to  be  considered  hetc  arc  the  materials  factor, 
mateiial  processing  factor  and  lubrication  laetoi  . 

All  of  the  rolling  traction  elements  in  the  Nasvytis 
■nuiitroiler  tiaeiion  diivc  arc  made  from  cousuiirable 
electrode  vaeuutn  melted  (CVM)  materials  (ref.  8). 
The  sun  roller  is  made  from  through-hardened  CVM 
A1SI-,‘>2100  bearing  steel  w  ith  a  Roekwcll-C  hardness 
of  61  to  63.  According  to  reference  5,  in  a  bearing 
application,  the  use  of  CVM  A1S1-521(K)  steel  would 
result  in  a  eoinbincd  material  aird  processing  life 
imtltiplii-aiion  factor  of  6.  The  other  ti action 
components  arc  mamifactuied  from  CVM  SAE-9310 
steel  that  is  case  carburized  to  a  RocKweli-C  batdncss 
of  60  to  63.  Reference  5  suggests  that  the  use  of  iliis 
material  will  give  a  life  multiplication  equivalent  to 
through-hardened  CVM  AlSl-52100. 

It  is  well  known  that  the  thickness  of  Inbricani  film 
which  separates  contacting  machine  elements  can 
have  a  strong  influence  on  the  contacts’  fatigue  life. 
In  reference  5,  the  effectiveness  of  a  lubricant  film  in 
terms  of  the  ratio  of  film  thickness  to  composiic 


siiilace  lougltness  /r/o  is  iclatcd  to  fatigue  life.  T  lie 
compitsitc  suifacc  toughness  is  defiticd  by 

•’  ^  -*  "Ji 

vvlrcic  arid  ('/)  ate  ilic  smfacc  finishes  of  the 
matitrg  bodies.  The  lubricant -fatigue  factoi  has  been 
assigned  a  value  of  1 .1)  whcti  /i/u  has  a  value  ol 
apmoximatcly  1.3.  Details  of  the  film  thickness 
calculation  using  the  Aichaid  and  Cowking  foinuila 
lecoinmcndcd  in,  rclcicncc  5  ate  given  in  icfcicncc 
11.  The  lubiicant  used  in  the  Nasvyti'.  nmltitollci 
traction  drive  in  reference  8  was  a  synthetic 
cycloaliphatic  traction  lubiicant,  and  the  roMci 
surface  roughnesses  were  less  than  0.2  nm  (8  /bn  ) 
rms.  A  contact  inlet  leinpcraiurc  of  344  k  (160^  1 ) 
was  assumed.  Due  to  inlet  shear  heating  and 
starvation  effects,  it  is  unlikely  that  the  minimum 
run  flint  thickness  would  exceed  the  composite 
surface  roughness  by  moic  titan  a  factor  of  4  (icf. 
16).  Accordingly,  the  li/o  latio  in  this  invcsiigatioii 
will  be  limited  to  a  value  of  four. 

One  additional  factor  not  treated  in  lefcience  5  but 
iiitporiaiti  to  iractioit  drive  contacts  is  the  potentially 
dcleteriou.s  effect  which  iiaction  itself  has  on  the 
lolliiig-elcmcni  fatigue  life  of  ilte  contacting 
clcmcitiv.  Some  iitvesiiguiions  (icf.  17)  !•■ave  found  a 
decrease  in  life  in  rolling-elciticiii  fatigue  tests  in 
wliich  some  sliding  or  iraciion  is  introduced.  l  aiigiic 
tests  (ref.  18)  vvii)i  increased  spin  (i.e.,  rotational 
sliding  within  the  com  act  area)  also  showed  a 
reduction  in  fatigue  life. The  analysis  of  reference  2(1 
has  shown  that  tangential  or  traciion  contact  forces 
reduce  the  depth  at  which  the  maxiiiium  shearing 
stress  occurs.  Since,  in  the  l.unabcrg-Palingrei. 
tlicory  (lef  10),  life  ilieoiciically  dccieascs  wiin  a 
decrease  in  the  depth  to  the  ntaximuii'  orthogonal 
reversing  slicai  stress  (sec  eq.  (3)),  traciion  forces 
may  .ict  to  reduce  contact  fatigue  life.  To  account  lot 
this  rediiciion  air  addiiional  life  adjustment  factor, 
the  Iiaction  faciin  ,  will  be  iniroduccd.  Given  the  fact 
that  little  conclusive  traction  contact  fatigue  life  data 
is  available,  this  facioi  for  the  purposes  of  this 
analysis  will  be  arbitrarily  given  a  value  Tot  all 
coiuacis  of  0,5.  This  factor  can  be  tieaied  as  an 
arbiiiary  bit  of  conservatism  unit!  mote  definiiive 
values  can  he  deici mined. 

A  suiiimary  of  the  life  adjustment  faetoi  s  appears 
in  table  1  for  ilii  ee  typical  points  which  cover  a  range 
of  operating  conditions  up  to  the  maximum  rated 
condition  of  194  kilowatts  (200  tip)  at  75  (KK)  rpni 
input  speed.  It  is  apparent  fioin  table  1  iliat  these 
factors  vary  relatively  little  ovei  the  entire  operating 
condition  spectruin.  Also  listed  iit  table  !  is  the 
l.undbeig-Palmgren  hfe  at  each  contact,  the  contact 
life  adjusted  for  life  factors  and  the  adjusted  drive 


Tabu  i  -si'NImarv^))  ct)NiA(.i  ionihiions.  i  ii i  adjusimini  i  ach)RS 

AND  DKIVl  l.ll  L  lORlHRll  DRIVI  Ol'l  RAl  INU  (.  UNUI I  IONS 


Sun- 

lirst- 

Second 

Sun 

1  list- 

Se».kmd 

Sli'i 

1  nsi 

Sectmd 

fiisi 

secs'tid 

planet- 

lust 

second 

planet 

til  St 

sectvnd 

planet  - 

planet 

planci 

rum 

planet 

planei 

firm 

planet 

planei 

inili: 

Ift.ti  kW  (22. 

:  lip) 

74. ft  kW  (UK)  hp) 

149  kW  (2WI  lip) 

25  000  r  pm  sun  speed 

75  IKK*  rpni  sun  speed 

75  OW  1  ptn  sui 

speed 

Maicnal  and  prc'iccbsiiiji  t'acioi 

6 

6 

6 

6 

fl 

f. 

Suifacc  sp«d.  111  see 

3(1  ft 

■  14  3 

14  5 

no 

43  0 

43  5 

no 

4V0 

4V5 

Normal  load,  N 

!3.M 

229ft 

1470 

2040 

3440 

2940 

40K0 

ft890 

Maxiiiiuni  coiiiao  pressure,  Cit’a 

1  11 

1  18 

0.950 

1  27 

1  35 

1  (>9 

1  ftO 

1  70 

IJ? 

I'iim  thickness,  inn 

1  14 

0  489 

0.972 

1  14 

1.07 

1.14 

1  14 

i.o: 

1  14 

h'O 
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2  ftS 

1  8S 
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2  fts 

2  ft5 
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2  h5 
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0  5 
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0  5 

0  5 

0  5 

0.5 

0  5 

0  5 

Total  lite  adju^imcni  factor 

7.95 

5.55 

7.65 
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7.80 
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7.95 

7.80 

7.95 

1  undliCrg-Palriiiireii  life,  tir 

Ii  4(Vt 

4(M() 

1.41  X  ll)0 

1130 

399 

1 39  (88) 

141 

40.9 

17  400 

Admsted  life,  lu 

6tXl 

22  4(8) 

1  08x  lO"^ 

8980 

3110 

l.)l  X  |06 

1120 

389 

138  (88) 

Adjusted  drive  life,  hr 

18  8(81 

2440 

305 

system  life.  Figure  9  shows  ilie  life  over  a  range  of 
ittpui  speeds  and  power  levels  based  on  the  contact 
fatigue  life  theory  and  life  adjustiueiit  factors. 

Effect  of  Relative  Rolici  Si/c  oii  Drive  Life 

The  effect  of  relative  roller  size  on  Nasvyiis 
traction  drive  system  life  can  be  studied  using  the  life 
analysis  described  previously.  For  a  given  ratio, 
number  of  planets  and  an  arbitrary  overall  size,  the 
relative  roller  diameters  can  be  varied  to  obtain  the 
maximum  total  drive  system  life. Once  the  cptiimim 
relative  roller  diameter  ratios  have  been  determined, 
one  can  scale  the  overall  drive  package  size  to  meet 
the  given  torque,  speed  and  service  life  requirements. 
The  procedure  to  be  followed  is  simply  to  select 
allowable  roller  diameters  and  center  locations  and 
then  calculate  the  latigue  life  for  an  arbitrary  liug 
roller  size.  The  roller  center  locations  arc  defined  by 
dimensions  B  and  C  of  figure  5,  and  the  following 
equations: 

B  =  r(,-rs 

C  =  zi+r2 

The  allowable  values  of  B  and  C  are  determined  from 
geometric  constraint  relations.  The  most  obvious 
constraint  is  that  the  drive  must  have  a  specific  ratio. 
The  first  row  of  rollers  must  be  of  a  diameter  that 
will  tit  without  mutual  interference  according  to  the 
selected  center  location  dimension  C.  The  same  is 
required  of  the  second  row  rollers  and  dimension  B. 
First  row  planet  (driving  track)  radius  rj  was 
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coiisi rained  to  be  no  Mnallci  than  5  percent  of  the 
ring  bore  radius  /■(,  for  bending  strength 
considerations  of  tire  first  row  roller  shaft.  In 
addition,  it  is  required  tliat  a  toggle  action  operate 
between  the  tollers  in  order  to  maintain  sufficient 
normal  load  at  each  traction  contact .  Toggle  action  is 
the  amplification  of  notmal  forces  acting  between  the 
first  aiid  second  rovt  rollers  when  the  second  row- 
rollers  wedge  themselves  between  the  first  row 
rollers.  This  toggle  action  or  wedging  effect  helps 
match  the  imposed  normal  load  to  that  required  to 
transmit  the  torque  for  a  given  traction  coelficient. 

Figure  10  shows  the  contraints  and  region  of 
allowable  roller  center  locations. In  order  to  eliminate 
drive  size  front  tlic  optimization,  the  axes  in  this 
figure  have  been  nondimensionalized  by  dividing 
dimensions  B  and  C  by  the  ring  radius  r^.  Each 
constraint  is  in  the  form  of  a  boundary  which 
separates  a  region  of  acceptable  (B/r^,  C/r()  values 
from  those  that  are  not.  In  figure  10,  the  values  of 
(B/zf,,  C/rf,)  which  meet  all  the  given  constraints  lie  in 
the  cross-hatched  region  as  shown.  A  drive 
configuration  having  values  of  (S/r^,,  OTf,)  which  lie 
outside  this  allowable  region  vould  have  one  or  more 
of  tire  following  deficiencies:  toller  to  roller 
interference,  no  toggle  action,  radius  ^3  too  small  or  a 
ratio  not  equal  to  tlte  chosen  ratio  of  15  to  1.  For 
each  allowable  combination  of  center  locations  C/>(t 
and  B/Zft,  there  is  a  range  of  allowable  planet  and  sun 


lijun-  10.  -  I’pt  inissiDli' region  for  roller  centej  loialion 
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roller  dianreters.  That  is,  the  diameter  of  the  sun  may 
be  made  larger  if  the  diameters  of  the  rcmaitiing 
rolleis  arc  adjusted  in  size  consistent  with  the  overall 
desired  15-to-l  ratio.  The  range  of  such  adjustments 
decreases  v/ith  increases  in  B/r(,  for  a  fixed  value  of 
C/ff,,  reaching  the  limiting  condition  at  the  right 
most  boundary  (ratio  constraint  boundary)  where 
ihcre  IS  only  one  set  of  toller  sizes  that  yield  the 
desired  ratio. 

The  calculations  showed  that  for  any  selected  point 
(B/r^,  C/rt)  the  life  was  greatest  if  the  maximum 
allowed  sun  diameter  was  used.  The  variation  in  life 
with  sun  roller  dlanieier  (with  fixed  B/Zf,,  C/r^)  was 
small;  typically  less  than  10  percent.  The  variation  in 
life  with  roller  center  locations  was  larger.  Figure  11 
shows  the  maximum  life  for  location  (B/Zf,,  C/Zf,) 
plotted  as  a  function  of  (B/Zf,)  and  (C/zf,).  The 
maximum  life  occurred  with  the  following 
nondimciisional  roller  diameters,  sun  ihrough  ring, 
respectively:  0.165,  0.234,  0.095,  0.326,  0.326,  and 
1.0.  The  maximum  life  occurred  at  (B/Zf,, 
C/z^)  =  (B/Z(„  C/Zf,)*  =(0.6733,  0.399)  which  iics  on 
the  ratio  constraint  boundary.  The  Nasvytis 
niuliiroller  traction  drive  reported  in  reference  S, 
which  had  a  slightly  lower  ratio  of  14.7  to  1,  was 
nearly  optimized  with  (B/Z(„  C/z^,)  =0.663,  0.387, 

The  life  vaiied  greatly  along  the  ratio  constraint 
boundary  line.  It  is  interesting  to  investigate  the 
reason  fci  the  decrease  in  life  as  the  geometry  is 
changed  from  the  optimum.  Figure  12  shows  the 
variation  in  drive  life  and  component  lives  with 
location  on  the  ratio  constraint  boundary.  It  should 
be  rectnphasized  that  all  roller  diameter 
combinations  on  this  figure  yield  a  constraint  ratio  of 
1 5  to  1 .  The  rnaxinium  life  is  apparently  achieved  as  a 
trade-off  -among  the  component  lines.  The  sun  atid 
driven  tract  of  the  first  row  rollers  increase  in  life 
while  the  remaining  components  decrease  in  life  as 
B/Z(,  and  increase. 

F'ig’urc  13  shows  the  variation  in  toniact  presstirc 
with  loeation  on  the  ratio  constraint  boundary.  The 
stress  continually  decreases  on  the  sun-planet 
contact.  The  oihei  two  contacts  see  similar  stress 
trends  with  fairly  constant  levels  of  stress  up  to  the 
larger  values  of  C/rf,)  where  the  stresses 

rapidly  increase.  As  a  measure  of  the  overall  severity 
of  stress  the  rescaled  product  of  o^^,,  Opp,  and  Op^  is 
shown.  The  minimum  of  this  function  occurs  ai 
B'zi,  =  0.6773  whereas  the  maximum  life  is  at 
B*/r,s  =  0.6733.  The  trends  of  contact  stress  product 
along  the  ratio  constrain  line  correlate  well  with  life 
trends.  The  overload  factor  as  a  function  of  position 
on  the  ratio  constraint  boundary  was  also  inspected. 
For  values  o.'  (B/z^,,  C/rf,)  less  than  (B/rt,,  C/r(,)‘  the 
overload  factor,  or  the  ratio  of  applied  normal  load 
to  required  normal  load,  ranged  from  1.0  to  1.25. 
Above  (B/z.,,  C/rpY  the  overload  factor  increased 
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fi'jure  12.  ■  V'jrijtion  in  Oriv.*  component  live*,  js  lunclion  of  position  on  rotio  conslroinl  liouncldry  (ratio  ^ 
b  to  1'. 

exponentially  for  the  planet/planet  and  the  by  loss  of  toggle  action  and  the  resultant  need  for 

planct/ring  contact.  This  is  because  of  loss  of  the  high  overload  factors  on  the  planet-planet  and 

toggle  action  or  wedging  effect  between  planet  planet-ring  contacts. 

rollers.  In  summary,  the  loss  of  life  for  the  smaller  During  the  course  of  this  optimization  study,  both 
values  of  (B/rfi,  C/r^)  was  due  to  high  stress  on  the  the  design  traction  coefficient,  m°.  and  the  rollers’ 

high  speed  sun  element.  T!ie  loss  of  life  at  the  higher  transverse  radii  were  heiU  constant.  With 

val-es  of  C/rj)  was  due  to  high  stress  caused  adjustments  in  relative  roller  size,  the  peak  available 
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traction  coefficient  of  the  lubricant  will  vary 
somewhat  and  may  necessitate  a  small  change  in  pc’. 
However,  variations  in  would  be  expected  to  have 
a  relatively  minor  effect  on  the  results  of  this  study. 
With  regard  to  the  roller’s  transverse  radii,  it  is  clear 
tliai  changes  ifi  ihc  iransvcrsc  gcoHiCiry,  such  ss 
roller  width  for  line  contact  or  crown  radius  for 
elliptical  contact,  can  cause  significant  variations  in 
contact  stresses  and  contact  life.  In  selecting  the 
proper  transverse  roller  radii,  other  considerations 
such  as  the  effect  on  the  available  traction 
coefficient,  spin  velocities,  and  efficiency  should  be 
addressed. 

Another  constraint  on  this  analysis  is  that  the  basic 
drive  configuration  of  two  rows  of  five  planet  rollers 
each  was  not  varied.  However,  depending  on  the 
requirements  of  size  and  ratio,  a  drive  with  some 
other  number  of  planet  rollers  may  exhibit  a  higher 
drive  system  life.  For  a  truly  optimi  ted  drive  design 
method,  the  traction  coefficienct  and  transverse 
contact  geometry  would  be  varied,  and  the  entire 
procedure  repeated  for  different  numbers  of  rollers 
and  planet  rows  until  a  maximum  life  was  found. 


Concluding  Remarks 

The  analysis  presented  deals  with  only  one  of  tne 
factors  that  should  be  considered  when  designing  a 
traction  drive.  Although  fatigue  life  is  important,  the 

......  _r  -_n _ 
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should  also  be  addressed.  Since  traction  drives, 
unlike  gearing,  rely  heavily  on  the  amount  of  traction 
that  can  be  generated  within  the  contact,  the 
lubricant,  operating  conditions,  contact  geometry, 
and  surface  finish  must  be  carefully  selected  on  the 
basis  of  both  performance  and  life.  Although 
traction  drives  relative  to  other  drive  mechanisms  are 
remarkably  simple  in  principle,  the  technology 
involved  in  arriving  at  the  best  design  is  rather 
sophisticated. 


Summary  of  Results 

A  contact  fatigue  life  analysis  for  a  Muliiroller 
Planetary  Traction  Drive  was  performed.  The  life 
analysis  takes  account  of  stress,  stressed  volume,  and 
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depth  of  occurrence  of  the  critical  stress.  The 
methodology  is  based  on  a  modification  of  the  well- 
known  Lundberg-Palmgren  theory.  A  design  traction 
coefficient  of  0.05  was  maintained.  Life  adjustment 
factors  were  included  in  the  analysis  based  on 
advancements  in  rolling-element  bearing  technology 
in  the  area  of  matciials,  processing,  lubrication  and 
design.  Elastohydrodynamic  film  thickness 
calculations  were  made  for  representative  operating 
conditions,  and  the  possible  adverse  effects  of 
traction  on  rolling  contact  life  were  considered. 

The  following  results  were  obtained  from  a  study 
of  a  14.7  to  1  ratio,  Nasvytis  Multiroller  Traction 
Drive  with  two  rows  of  five-stepped  planet  rollers 
each  in  a  9-kilogram  cluster  of  approximately  0.21 
meter  diameter  and  0.06  meter  width  having  sun 
roller-first  row  roller  (driven  track)  contact,  first  row 
roller  (driving  track)-second  row  roller  (driven  track) 
contact  and  second  row  (driving  irack)-ring  roller 
contact  ratios  of  1.28,  3.87,  and  2.97,  respectively: 

1.  The  total  drive  system  Lundberg-Palmgren 
10-percent  fatigue  life  with  life  adjustment  factors 
ranged  from  18  800  hours  at  16.6  kilowatts  (22,2  hp) 
and  25  000  rpm  sun  roller  speed  to  305  hours  at  149 
kilowatts  (200  hp)  and  75  0()0  rpm,  with  a  life  rating 
of  2440  hours  at  a  nominal  operating  condition  of 
74.6  kilowatts  (iOO  hp)  and  75  i'kXi  rpm. 

2.  For  a  given  ratio  and  number  of  planets  an 
optimum  combination  of  roller  diameters  and  center 
locations  exists  for  maximum  drive  system  life. 

3.  The  shortest  lived  contact  was  the  first  planet- 
second  planet  contact  since  the  highe.si  maximum 
Hertz  stress  occurred  there. 

4.  For  a  given  fatigue  life,  traction  coefficient  and 
drive  configuration,  torque  capacity  is  proportional 
to  the  2.8  power  of  overall  drive  size.  For  elliptical 
contact  shapes  and  constant  torque  drive,  life  is 
proportional  to  the  8.4  power  of  size. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  August  29,  1980. 
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